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Abstract 

We analyze the coupling of CP-even and CP-odd Higgs bosons to a photon and 
a Z boson in extensions of the Standard Model. In particular, we study in detail the 
effect of charged Higgs bosons in two-Higgs doublet models, and the contribution 
of SUSY particle loops in the minimal supersymmetric extension of the Standard 
Model. The Higgs-yZ coupling can be measured in the decay Z —> 7-I-Higgs at e+e“ 
colliders running on the Z resonance, or in the reverse process Higgs ^ Z 7 with the 
Higgs boson produced at LHC. We show that a measurement of this coupling with 
a precision at the percent level, which could be the case at future e^e“ colliders, 
would allow to distinguish between the lightest SUSY and standard Higgs bosons 
in large areas of the parameter space. 


1. Introduction 

The study of the electroweak symmetry breaking mechanism is one of the most im¬ 
portant goals of present and future high-energy colliders. Once the hrst evidence for 
Higgs particles is established, it will be crucial to ascertain the underlying dynamics of 
the Higgs sector. This can be achieved by measuring the couplings of the Higgs particles 
to the other fundamental particles: in the Standard Model (SM), fermions and gauge 
bosons acquire masses through the interaction with the Higgs field and the size of their 
couplings are set by the masses. This is a fundamental prediction of the Higgs mechanism 
which has to be tested experimentally. 

Among these Higgs couplings, the couplings to photons are interesting in many re¬ 
spects. First, the interaction of the Higgs particle with photons does not occur at the 
tree-level since the photon is massless. The Higgs-photon coupling is therefore induced by 
loops of heavy charged particles. In the SM, this occurs via W boson and heavy fermion 
triangle loops. Since the couplings of these particles to the Higgs boson grow with the 
mass, they balance the decrease of the triangle amplitude with increasing loop mass, and 
the particles do not decouple even for very large masses. Therefore these processes can 
serve to count the number of heavy charged particles which couple to the Higgs boson. 

In supersymmetric (SUSY) theories, the Higgs sector must be extended to contain (at 
least) two doublets of scalar fields, leading to the existence of five Higgs particles: two 
CP-even Higgs bosons h and FT, a CP-odd Higgs boson A as well as two charged Higgs 
particles |^. Depending on the value of tan (3 [the ratio of the two vacuum expectation 
values of the Higgs fields], the lightest SUSY Higgs boson h is constrained to be lighter 
than Mh ^ 80-130 GeV in the minimal version. For small values tan/5 ~ 1.5, which 
are favored by Yukawa coupling unification [Q, the Higgs boson has a mass which does 
not exceed Mh ~ 80 GeV, and therefore can be produced at LEP2 |^. In the decoupling 
regime Ma Mh ~ Mh± > Mz i, the h boson has practically the same properties as 
the SM Higgs particle; the MSSM and SM Higgs sectors then look almost the same, and are 
very difficult to be distinguished. However, additional contributions to the Higgs-photon 
couplings will be induced by loops with charged Higgs bosons, charginos and sfermions. 
Since the SUSY particles do not couple to the Higgs boson proportionally to their masses, 
their contributions decouple for high masses. If, however, some of these particles are not 
too heavy their contributions can be large enough to allow for a discrimination between 
the lightest SUSY and standard Higgs bosons even in the decoupling regime. 

The Higgs coupling to two photons has received much attention in the literature 
[^, 1^. In the SM, the coupling of the Higgs boson to a photon and a Z boson 
has been discussed in Ref. [§. If Mz < Mho < 130 GeV, the H^Z'y vertex can be 
measured in the decay process ^ Zj. At LHG the production rate for light Higgs 
bosons is very large, a{gg —> H^) ~ 100 pb [Q, and despite of the small branching ratio 
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Z'y) ~ 10“^, one would still have 0(10^) Z'j events if the luminosity is high 
enough, C ~ 10^‘^cm“^s“^. The coupling could be therefore measured if the background 
can be reduced to a manageable level and if the theoretical prediction for the cross section 
is well under control [^. 


If Mjjo < Mz, the H^Z'j coupling can be measured in the reverse decay Z H^'j. 
At LEP and SLC, the rates are however rather small, BR(Z —< 10“^, leading to 
only a few events. However, at future e’''e“ colliders 0 with the expected integrated 
luminosities of / £ ~ 50 fb“^, running a few weeks on the Z resonance would allow to 
obtain a very large sample of Z ^ events. A very precise measurement of the H^Zj 
coupling would be possible in this case. If the Higgs boson is discovered at LEP2, one 
would then use the Next Linear Collider to measure the H^Zj coupling and check 
whether the Higgs boson is SM-like or not. This measurement would be then equivalent 


to measuring the coupling at high-energy yy colliders 


In supersymmetric theories, the couplings of the light CP-even and CP-odd Higgs 
bosons to Zy have been studied some time ago in Ref. |jT^ (see also fl^)- However, 
these analyses need to be updated for several reasons: (i) the radiative corrections in the 
MSSM Higgs sector turned out to be very large @], and therefore must be included; (ii) 
the Higgs couplings to top squarks can be strongly enhanced if squark mixing is included 
and this might induce large contributions to the Higgs-Zy coupling, a possibility which 
has been overlooked; (hi) stronger experimental bounds on the masses of charginos and 
sfermions are now available []^, eliminating a large part of the SUSY parameter space 
where contributions from these particles are large; (iv) hnally, a fully analytic expression 
for the contributions of charginos and top squarks with different masses is still lacking. 


In this paper, we address all the previous points. Furthermore, we discuss in some 
details the possibility of using the Higgs-Zy couplings to discriminate between the Stan¬ 
dard Model and its extensions. In particular, we analyze to what extent one can use 
the SUSY loop contributions to distinguish between the standard and light SUSY Higgs 
boson in the decoupling limit where all other Higgs [and SUSY] particles are too heavy 
to be produced directly; we will show that this is indeed possible if the coupling can be 
measured at the percent level. We also show that in a general two Higgs-doublet model, 
the contributions of charged Higgs bosons do not necessarily decouple from the Higgs-Zy 
amplitude for large masses, contrary to the SUSY case. 


The paper is organized as follows. In the next section, for completeness and to set 
up the notation, we discuss the iL°Zy coupling in the Standard Model. In section 3, 
we analyse the Higgs-Zy coupling in the two Higgs-doublet extension of the model. In 
section 4, we discuss the various loop contributions of the SUSY particles to the coupling 
in the MSSM, paying special attention to the small tan/3 region and the decoupling limit. 
Our conclusions are given in section 5. 
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2. The Higgs—Z—photon coupling in the SM 

For the sake of completeness and to fix our notation, we first discuss the Higgs-Z 7 coupling 
in the Standard Model. The H^Z'y vertex is mediated by W boson and heavy quark [in 
practice only top and bottom quark] loops; Fig. la. It can be decomposed into the 
following tensorial structure: 

V[Z>^{p^),Y{P2),H\p,)] = Fopi^p’i + + Fgp^p^ + F,p^^p\ + Tge'^-^^pi^pa^ (1) 


For on-shell particles, only the form factor Fq contributes to the decay widths; normalized 
to e^/ (swMw) with = 1 — = sin^ 6w, it is given by 


Fq — Aw + Af = M| 


cot 9wFw + 2Q f Nc 

f 


m) li-2s^^Qf 
M| sw cw 



( 2 ) 


with Qf, Jg and m/ the charge, the weak isospin and the mass of the fermion /; iVc = 1 
for leptons and = 3 for quarks. 


In the following we discuss the various contributions in the case of the decay Z —> 
the amplitudes for the reverse decay can be simply obtained by crossing. In terms of the 
Passarino-Veltman three-point scalar functions 


Co, 2 {m‘^) = Co, 2 (M|, 0, m, m) 


( 3 ) 


the fermionic and W contributions^ are found to be 



Co(^/) + 4C2(^/) 


M^o 


(1 - 24) + 2(1 - 64) 


C'2(M^) + 4(1 - 44)C'o(M^) (4) 


Cq is the scalar integral, and C 2 is a short-hand notation for C 2 = Cii + C 23 where the 
expressions of the 4 can be found in [^. Since there is only one mass running in the 
loops, the functions C'o(^^) and C 2 {'ni?) have a rather simple form; in terms of the scaled 
variables tz = /M'^ and th = 4m^/M|fo, they are given by the known expressions [|l| 


4m4o(m2) = [f{Tz) - /(r^)] (5) 

Tz — Th 

with the functions / and g dehned by Jl| 


/4 = 


arcsm 


^14 


it: 


r > 1 
r < 1 


( 6 ) 


^We have calculated the amplitudes in the Feynman gauge; however, the results for the fermion and 
IF loops are separately gauge invariant if all external particles are on-shell. 
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9{'r) = 


\/t — 1 arcsinyi/r 


ITT 


r > 1 
r < 1 


(7) 


logT^ 

The W boson and top quark form factors and Af are shown in Fig. 2a as a function 
of the Higgs boson mass. The b quark contribution as well as the contributions of the 
other fermions are much smaller, due to the small masses. In the range of interest 70 GeV 
< Mffo < 130 GeV, the W contribution is by far dominant, being one order of magnitude 
larger than the top quark contribution; the two amplitudes interfere destructively. The 
QGD corrections to the top quark loop are small, being of 0{as/'n') . 

The decay rate for the process Z —>■ reads in terms of Fq: 




1927r4 


^ Ml 


1 - 




( 8 ) 


This rate, normalized to the total decay width Tz — 2.5 GeV, is displayed in Fig. 2b. The 
branching ratio varies from ~ 10“® for masses Mfjo ~ 50 GeV [which are ruled out in the 
SM, but are still possible [|T^] in extensions of the model] to ~ 10“^ for Mho ~ 80 GeV 
which can be probed at LEP2 [Q. This means that only a few events can be produced 
at LEPl with the present sample of (P(IO^) Z bosons collected by all four collaborations. 
However, a future collider with the expected yearly integrated luminosity of / £ ~ 100 
fb“^, will be able to produce (9(10^*’) Z bosons per year; this translates into (P(IO^) 
Z —> H^'y events for Higgs boson masses not too close to the Mho ~ Mz threshold. Since 
the signal is very clean [the photon being monochromatic and the decay products of the 
Higgs boson, bb, being efficiently tagged with micro-vertex detectors], one could 

measure the H^Z'y coupling with a statistical precision of a few percent allowing for a 
stringent test of the H^Z'y coupling. One would therefore check whether the coupling is 
SM-like, and measure with a good precision the H^WW and H^tt coupling. 

If Mho > Mz, the decay rate for the reverse process —>■ Zj reads 


F(77° 




647r"^ 


HO 


Ml ' 

Mi., 


(9) 


with Fq given by eq. (2). The branching ratio when the Z boson is decaying into electron 
and muon pairs is also shown in Fig. 2a [with BR(Z —>■ ~ 6%; other decays 

of the Z bosons will be rather difficult to extract from the background at the LHG]. In 
the mass range Mho ~ 120 GeV, the branching ratio is of the order of 10“^. With the 
(P(10®) Higgs bosons produced at LHG in the main production mechanism gg 
with an expected yearly luminosity of / £ ~ 100 fb“^, a few hundred —>■ Zj events 
could be collected in a few years of running, if background problems can be reduced to 
a manageable level. This rises the hope to measure the H^Zj coupling once the Higgs 
boson is observed in the —>■ 77 mode for instance. 
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3. The coupling in Two—Higgs Doublet Models 

In a Two Higgs-Doublet Model (THDM), the most general Higgs potential compatible 
with gauge invariance, the correct breaking of the SU(2)xU(l) symmetry and CP con¬ 
servation is given by |]l| 

V = Ai(|0ip — -|-A2(|02p —~ + (|02p — 1 ^ 2 )]^ 

+ A4[|0ip|02d — 10102^] + A5[Re(0402) — '^ 1 ^ 2 ]^ + Ae [Ini(0|02)]^ (10) 


with 01, 02 the two Higgs-doublet helds and Vi,V 2 their vacuum expectation values. 
We have also assumed that the discrete symmetry 0i —>■ —0i is only broken softly; an 
additional term, A 7 [Re( 0 i 02 ) — uiU 2 ]Ini( 0 i 02 ), can be eliminated by redehning the phases 
of the scalar helds |l|. Parameterizing the Higgs doublets by 


\Vi + T/i+iXlJ \V2 +V2 +iX2j 

one obtains for the mass terms in the CP-even Higgs sector 

/ 4(Ai Aa)^! -1- A 5 V 2 (4A3 -h A5)nir;2 
( 4 A 3 + A5)r;ir;2 4(A2 + A3)t;2' + 

while in the CP-odd and charged Higgs sectors, one has 



•^6(Xl) X 2 ) 


-V 1 V 2 

vf 



, A4(0i , 02 


-^1^2 


-^1^2 Vl j \X2/ \ -V 1 V 2 

Diagonalizing the mass matrices, one obtains the physical masses 

-|- A^22 i (-Adll — A422)^ + 4W1i2 
M\ = Aqv'^ and = A 4 n^ 



Mi. = 2 


( 11 ) 


( 12 ) 


(13) 


(14) 


with = vf + V 2 and Ai the mass matrix of eq. (12). The mixing angle a in the CP-even 
Higgs sector is obtained from 

Af 11 — Af22 . ... 2A4i2 


cos 2a = 


, sin 2a = 


\/ (Af 11 — M.22)‘^ + 4Af 12 (Af 11 — M.22)‘^ + 4AIi2 

Inverting these relations, one obtains the A’s in terms of the Higgs masses, and a,0: 

1 . 2 ^ ^2 2 ^ ^2-, sin 2 a 1 X 1 % — Ml A5 , sin^d 

[cos^ aMjj + sin^ aMl) ^ h , 5,, 


(15) 


Ai 

A 2 

A 3 


4 cos^ I3v‘^ 

1 

4 sin^ f3v'^ 
sin 2a Mjj 


sin 20 


-R« + h(i 

4ti2 4 ' 


sin^ aMjj + cos^ OiMl) 


sin 2 a Mjj 


cos^ 0 

,2 


), 


sin 20 


Ml ^ . _ cos^0 . 
4r;2 4 sin^ 0 ’ 


Ml 


sin 20 4n^ 


As 

T 




— —s- 


(16) 
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As one can see, the parameter A 5 can not be fixed by the masses and the mixing angles, 
unless one imposes a strict —>• —0i symmetry resulting in A 5 = 0, or by using the SUSY 
relation A 5 = Ag = as will be discussed later. 

In a general THDM, the four masses Mh-, Ma and Mh± as well as the mixing 
angles a and (3 are free parameters. The interaction of the Higgs bosons with fermions 
are model-dependent; here, we will consider the model where one Higgs doublet couples 
only to up-type quarks, while the other doublet couples only to down-type quarks and 
charged leptons [the so-called Model H [^] which occurs in SUSY models for instance]. In 
this case, the couplings of the neutral Higgs boson, collectively denoted by <h, to fermions 
and massive gauge bosons are given in Tab. 1, normalized to the SM Higgs couplings. 
Due to CP invariance, the pseudoscalar A does not couple to W and Z bosons. 


$ 

Q^UU 

g^dd 

g^vv 

h 

cos a/ sin /3 

— sin a/ cos [3 

sin(/? — a) 

H 

sin a/ sin 

cos a / cos [3 

cos(/? — a) 

A 

1 / tan/? 

tan/? 

0 


Tab. 1 Higgs boson couplings to fermions and gauge bosons in the THDM compared to 
the SM Higgs couplings. 


For the CP-even Higgs bosons, the couplings to photons and Z bosons receive contri¬ 
butions from W and top/bottom quark loops as well as contributions from charged Higgs 
boson loops (Fig. lb). The structure of the vertex is again given by eq. (1), and only the 
form factor Fq contributes for the decay. It is given by [4> = h,H] 


Fn = Mi 


cot 9w SrpVvFw + cot 9w g(l>H+H-FH± 


f 


mi J, 


2^2 


Qf 


^Ml 


SwCw 


Stilus 


(17) 


The functions Fw and Ff are the same as previously, while the function Fh± for the 
charged Higgs contribution reads in term of the C 2 function defined previously (see also 

H): 

Fh* = 4C2(MJ±). (18) 

The couplings g(f>vv and g^ff can be taken from Tab. 1, while the couplings of the CP- 
even neutral Higgs bosons to charged Higgs bosons in the THDM [using a normalization 
similar to the one for the W boson and the fermions] are found to be 

— Asu^ cos(/3 -I- a) 


ghH+H- — 


gHH+H- — 


Ml- 


Ml 


M'i, 


sin 2,3 
Asu^ sin(/3 -|- a) 
sin 2(3 


+ 


2M^± 


Ml 


+ 


2M^ 

2Mjj± - Ml 


sin(/? — a) 


H 


2M2, 


cos(/d — a) 


(19) 
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In the limit of very heavy bosons, the C 2 {M‘jj±) function reduces to 

C 2 (Mi±) ^ 1/(24M|±) , 

while the coupling of the h boson [that we assume to be the lighter CP-even Higgs state] 
to the charged Higgs boson approaches the limit 

ghH+H- sm(/3 - a) 

assuming that ^'h±- This leads to a hnal contribution which is proportional 

to sin^(/5 — a) = g^vv- Therefore, the charged Higgs contribution to the hZy coupling 
in a general THDM does not decouple, contrary to the case of SUSY models as will be 
discussed later. However, the contribution is suppressed by the large factor 1/24 
and compared to the W boson loop [which in a THDM, is also damped by the factor 
sin(/5 — a) compared to the SM case], it is two orders of magnitude smaller. The decay 
widths r(Z —hy) or r(h — Zj), given by eqs. (8-9) with M^o will therefore 

be hardly sensitive to this loop effect. 

We now turn to the case of the pseudoscalar Higgs boson A. Due to CP-invariance, 
the AZ'j coupling is induced only by fermionic loops, since A does not couple to W and 
bosons. Its tensorial structure is given by the same expression as in eq. (1), but here 
only the form factor contributes. The decay widths r{Z ^ A'j) or r(H — Z'j) are 
given by eqs. (8-9) with M^o Ma and Fq replaced by 

2 

F, = -y 24,0;)g.,„C„(,q), (20) 

/ 

In the general THDM, a numerical analysis is rather complicated [and not very telling] 
since besides the four masses Mh, Mh-, Ma and Mh±, we have the mixing angles a and 
(3 as additional parameters, not to mention the parameter A 5 which is also independent. 
The Higgs-Zy couplings can vary widely compared to the SM coupling, although the 
dominant W boson amplitude is always suppressed by the factors sin(/5 —a) or cos(/5 —a) 
in the case of the CP-even bosons or absent in the case of the pseudoscalar A. To simplify 
the discussion, we will use the constraints provided by supersymmetry: in the MSSM, the 
Higgs sector is described at the tree-level only by two free parameters that we chose to 
be tan/3 and the pseudoscalar mass Ma- The masses of the CP-even Higgs bosons are 
given by 

Ml„ = i [m 5 + Ml T ] ( 21 ) 

while the mass of the charged Higgs boson is simply given by 

Ml^ = Ml+ M^ 
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( 22 ) 



The mixing angle a is related to Ma and tan (3 by 

tan2a = tan, -| < a < 0 (23) 

However, these relations are affected by large radiative corrections 0] which must be taken 
into account. We will therefore include the leading radiative correction to the Higgs masses 
and couplings which grows as and logarithmically with the common squark mass that 
we fix to 1 TeV, unless otherwise stated. In the MSSM, one has 1 < tan/3 < rrit/mi, from 
GUT restrictions, with the lower [tan/3 ~ 1.6] and the upper [tan/3 ~ 50] ranges favored 
by Yukawa coupling unihcation |^. We will mainly focus on the lightest CP-even Higgs 
boson, for which the maximum allowed value of the mass is about ~ 80 GeV for 

tan/3 ~ 1.6, and the particle is therefore accessible in Z decays. In the high tan/3 range, 
the maximal h mass can reach values ~ 130 GeV, and the decay h Zy would be 
kinematically possible. 

The W boson amplitude Aw = cot 6w QhvvFw is shown in Fig. 3a as a function of 
Mh for the three values tan/3 = 1.6,5 and 50. For low h masses. Aw is suppressed 
compared to the SM value, the suppression being more effective with increasing tan/3; 
in fact, for tan/3 ~ 50 the W contribution almost vanishes. With increasing h mass. 
Aw approaches the SM value which is reached for Mh = The sum of the top 

and bottom loop contributions is also displayed in Fig. 3a. Except when Mh ~ 
where the the form factor Af becomes SM-like, the t contribution is suppressed by a 
factor ghtt 1/tan/3, while the b contribution is enhanced by the factor ghbb ~ tan/3. 
Therefore, the t contribution is dominant for low tan /3, while for large tan /3 values the b 
contribution [which has opposite sign compared to the SM case] is strongly enhanced and 
becomes dominant. 


The contribution of the charged Higgs boson loop Ah± = cot 6w ghH+H-Fn^ is shown 
in Fig. 3b as a function of Mh± for the values tan/3 = 1.6,5 and 50. The behavior can 
be understood by recalling the expression of the ghH+H- coupling in the MSSM 


QhH+H- 


sin(/3 — a) + 


cos 2/3 sin(/3 + a) 


2^2 


(24) 


in which the radiative correction must also be included; see for instance Ref. [0 . For small 
Mfh± implying small Mh, the coupling is strongly suppressed for large tan/3 values and 
the contribution Ah± is small. For tan/3 ~ 1, the suppression is rather mild and the 
contribution can be large, reaching a few percent of the W contribution for Mh± ~ 100 
GeV. Gontrary to the THDM, Ah± decreases with increasing Mh± since in the MSSM, 
the ghH+H- coupling does not scale like the charged Higgs mass, and the contribution is 
damped by a factor 1/M^±. The charged Higgs boson therefore yields small contributions 
to the hZ'y coupling and decouples from the vertex for high masses. 
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4. The Higgs—Z—Photon coupling in the MSSM 

In the Minimal Supersymmetric extension of the Standard Model, the couplings of the 
CP-even Higgs bosons to the photon and the Z boson receive, as in the Two Higgs- 
Doublet Model, contributions from W bosons, top+bottom quarks and charged Higgs 
bosons. Extra contributions also come from charged supersymmetric particles: sleptons, 
squarks and charginos (Fig. Ic). The decay widths are again given by eqs. (8-9) with 
Mjyo —>• M^, and the form factor Fq reads: 


Fo 


Ml 


cot 9w QcjiVvFw + X! 2^/ — - wRl Pj, 

f tWz SwCw 

+ cot 6w g,f>H+H-FH± + cot 9wF^+ + ^ NcQjFj . 

f 


(25) 


The amplitudes from the fermions, W and bosons are, as in the previous section, 
given by eqs. (4) and (18), while the chargino contribution reads 


m. 


^x+ = H -Tr~ f ("^x+ ’ ^xt ’ ^xt 
j,kF {,2 Mw ^ 


m,n=L^R 

The couplings of charginos to the Z bosons are given by 

1 


) i: s™ al*,- 


^ZXjXk ^<f>xtXj 


(26) 


sh*,,- = - = - (^“^t + - Siisl) (27) 

while the couplings to the higgs bosons read 


9(j,y+y- Q jF<t> ^ji^4> 1 Ssy+Y- Sijd, 


FxJ x] 


<t>Xj Xj 


(28) 


with Ch/dh = sin a/ cos a and cn/dn = — cos a/ sin a. The elements Qij/Sij, as well as 
the matrices V and U which diagonalize the chargino mass matrix can be found in p0|. 

The function / entering the chargino form-factor is given by 


/(mi, m 2 , m 2 ) 


-2 


Co{mi,m 2 ,m 2 ) + (Pi (mi, m 2 , m 2 ) + 2(5*2 (m-i, m 2 , m 2 ) 


+2C2{m2,mi,mi) - (Pi(m 2 , mi, mi) 


(29) 


where (Pi and C 2 have now a more complicated structure since there are two particles 
with different masses inside the loop. In terms of the scalar Passarino-Veltman functions 
Hq, Bq and (Pq one has 


(Pi (mi, m 2 , m 2 ) 


(Pii(M|, 0, Ml] mi, m 2 , m 2 ) 

Bo{Ml] mi, m2) - Bo{Ml] mi, m2) 
Ml - Ml 


(Po(M|, 0, Ml; mi, m2, m2) 
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C2{mi,m2,m2) = 


Ci2(M|, 0, Ml] mi, m 2 , m 2 ) + C23(M|, 0, mi, m 2 , m 2 ) 

^{M^- Afffi "^1’ "^ 2 ) - ^o(M|; mi, m 2 ) 

+7^7772—-772V772 C'o(M|, 0, M^; mi, m 2 , m 2 ) 


' 2{Ml - Ml)Ml[ ^ 

+ (m2 - mi)5o(M^;mi,m2) + Ao(mi) - Ao(m2) • (30) 

The expressions of the scalar one-, two- and three- point functions Aq, Bq and Co are 


ylo(m) = m^ 1 — log 


Bo{p^,mi,m2) = 


, mim 2 m? — m^ , m 2 
2 - log-^ ^-log — 

mi 

^ ^2 _|_ J^2 _ p2 _|_ ^l/2^p2^ ^2^ ^^2^ 

2mim2 


Co{Ml, 0, Ml, mi, m2, m 2 ) = 
Ml - M| ^ 


ml — mf + + (tAi/2 (M?, mf, ml 


p is the renormalization scale, and the ultraviolet poles in Aq and Bo are subtracted since 
the amplitudes are finite; A is the usual two-body phase space function: \{x,y,z) = 
— 2{xy + xz + yz). 

Finally, the contribution of the squark and slepton loops to the hZ'j couplings reads 

Ff- = -8 E (34) 

j,k=l,2 

with the function C 2 defined in eq. (30). The squark couplings to the Z boson, including 
mixing between left- and right-handed sfermions, are given by 


9zfJi 


9zf2f2 


9zhf2 


SwCw 

1 

SwCw 


SwCw 


{li - Q/S^) cos^ Of- Q/S^sin2 Of 
-Qf sly cos^ Of + {li - Qfsly) sin2 Of 


sin 9f cos9f 


The mixing is proportional to the fermion mass, and in practice is non-negligible only for 
the partners of the third generation fermions. 
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The couplings of the Higgs bosons to squarks have a more complicated structure 
because of the squark mixing. In the case of the light CP-even Higgs boson h, they read 


9hhh = CIl cos^ Of + sin^ Of + 20^^ cos Of sin Of 
9hhh = C'Lcos2 0/ + C''*^sin2 0/-2C'^iCos0/sin0/ 

9hhh = C'Ri(cos2 0/-sin2 0j) + (C'^^-C']^^)cos0/sin0/ (36) 

with 


CL = 

f 2 N 

\h QfSw)9hvv M‘2 9hff 


^RR — 

{QfSw)9hVV J ^2 9hff 
z 


^RL — ~ 

2 M2 MHff] 

(37) 


where H/ is the soft-SUSY breaking trilinear term and /i the Higgs-higgsino mass pa¬ 
rameter; the couplings gtff and gcpvv are given in Tab. 1. The heavy CP-even Higgs 
boson couplings to squarks can be obtained from the previous ones, by performing the 
substitutions gh... ^ gH...i sin a —cos a and cos a —sin a. Note that in the case of the 
partners of the light fermions, the mixing angles and the fermion masses can be set to 
zero, and these couplings simplify considerably. 

For the pseudoscalar Higgs boson H, only top-|-bottom quarks and the charginos are 
contributing to the AZ'-j amplitude because of CP-invariance. The form factor has to 
be replaced by 


Fs = -M| 2gj Ye ^ - cot 9wF^, 

I j IVIz SwCw 

with the chargino contribution 
m. 

’^3 ( 

— ni m. 4.. 17^ , . .,. , 


Fx* = T. {s^}x; + Atx-)i!>Axtx; -<tx; 

j,k=l,2 ^ 


(38) 


(39) 


The gAx+x~ conplings are given by 


dx*x; “ ^ ~ ^ ’ ^*xtx- “ ^ F S,j cos P , 


(40) 


and the new fnnction g reads 


g{mi,m2,m2) 


-2 


Co(mi,m2,m2) Ci(mi,m2,m2) -h Ci(m2,mi,mi) 


(41) 


In the MSSM, the CP-even Higgs boson H is always heavier than Mz and therefore 
only the decay H Z'j is kinematically possible. However, even before allowing the 
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Z boson to decay into charged leptons, the branching ratio BR(if — 7 Z) is very small 
and the process will be very difficult if not impossible to be seen at the LHC. This is 
also the case of the decay A Zj. The process Z Aj will be possible if Ma ^ 80 
GeV, but then the pseudoscalar Higgs boson can be discovered at LEP2 in the associated 
production mechanism > hA and its properties can be studied. In the numerical 

analysis, we therefore focus on the light CP-even Higgs boson h and study in particular 
the low tan f3 scenario in which < 80 GeV and the h boson can be produced at LEP2 
in the process e^e“ — hZ. We will pay special attention to the decoupling limit where 
the h boson mimics the SM Higgs particle and the measurement of the hZ'j coupling at 
future e’''e“ linear colliders running at the Z resonance could help discriminating between 
the SM and MSSM scenarios. 


The contributions of the W, SM fermion and charged Higgs boson loops to the hZy 
have already been discussed in the THDM with the MSSM constraints. As discussed 
previously, for large values of Ma, the W and quark contributions are as in the SM, while 
the charged Higgs boson decouples and its contribution is negligible. The contributions 
of the slepton and the scalar partners of the light quarks, neglecting sfermion mixing, 
are shown in Fig. 4 as functions of the masses and for the three values tan/3 = 1.6, 5 
and 50 with Ma hxed to 1 TeV. We have summed over all slepton and squark [except 
stop] contributions, and used common masses mi and mq. As in the case of the charged 
Higgs boson, slepton and squark loop contributions to the h Z 7 decay width are very 
small, except when these particles are very close to their allowed mass values |^. For 


loop masses above 150 GeV, the current experimental bound on squark masses, they do 
not exceed the level of a few permille of the dominant W contribution and will therefore 
hardly be detected. 


The contribution of the top squark loops to the hyZ vertex depend on the soft SUSY- 
breaking scalar masses and mi^ [which are taken in general to be equal], as well as 
on the soft-SUSY breaking trilinear term A^, the Higgs-higgsino mass parameter p and 
tan/3. These parameters also determine the masses and mixing angle of the scalar top 
quarks and their couplings to the Higgs bosons. In the decoupling limit, Ma 3> Mz the 
contribution only depends on the combination = At — fi/ tan/3 and which can 
be traded against the mass of the lightest top squark In Fig. 5a, we show contour 
plots in the (jrit^, mij plane for which the contribution Ai, which includes the amplitudes 
of both top squarks, is 0.5, 0.2 and 0.1. For large Ma, i.e. in the decoupling limit, the 
amplitude Ai is symmetric for positive and negative mt^ values. For large \mt^\, the 
contributions are large and negative; for light enough top squarks, mi^ ~ 150 GeV, they 
can reach the value Ai ~ —.5 for |m^^| ~ 1 TeV, i.e. at the level of the top quark 
contribution. For a given , Ai is larger for higher values of because in this case 
the coupling htt ~ is strongly enhanced. For large mi^, the two top squarks will 
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have comparable masses and their amplitudes will partly cancel each other, leading to the 
contour = 0. For small there is a region around = 0 where no solution 

for mi^ < rrit is allowed when diagonalizing the mass matrix and this region is already 
excluded by CDF/DO data since there one has rriq < 150 GeV 0- The amplitudes in 
this region are positive and can reach values Ai = 0.2 which decrease with increasing top 
squark mass as expected. 

Fig. 5b shows the contribution Ai for = 100 GeV, i.e. away from the decoupling 
limit, and tan/? = 1.6 [solid lines] and 50 [dashed lines]. For low tan/? values, the 
symmetry around = 0 is lost; the picture is the same as in Fig. 5a for positive 
values, but the contribution A^ becomes smaller for negative For the high tan/? 

scenario, the contribution A^ becomes very small. 

The contribution of the charginos to the hZ'j coupling depends on tan /?, the gaugino 
mass parameter M 2 and the Higgs-higgsino mass parameter p [these parameters also £x 
the chargino masses]. The form factor A^ is shown in Fig. 6 in the {M 2 , fi) plane for 
tan/? = 1.6 and Ma = 1 TeV and Ma = 80 GeV. Gontours for \A^\ = 2,1,0.5 and 
0.2 as well as the region of the parameter space for which the lightest chargino mass is 
larger than 70 GeV [which approximately corresponds to the current experimental limit] 
have been drawn. The chargino contributions are rather large close to the M^+ = 70 
GeV boundary, reaching values A^+ ~ 1, and become smaller when one moves away from 
this boundary. However, in a large part of the {M 2 , parameter space, the chargino 
contribution is larger than = 0.2 and does not strongly depend on Ma- In fact, in 
this case, A^+ > 0.1 in the entire parameter space M 2 , n < 500 GeV which leads to a 
deviation of the h Z'y coupling by more than one percent. 

5. Summary 

We have analyzed the Higgs-Zy coupling in the Minimal Supersymmetric extension of 
the Standard Model. We have included the large radiative corrections in the Higgs sector, 
updated the contributions of the charginos and the top squarks to the coupling, and given 
fully analytic expressions for these contributions. In Two Higgs-Doublet Models, we have 
shown that the additional contribution from charged Higgs bosons do not necessarily 
decouple from the amplitude, if the mass is large. 

We have payed special attention to the case of the MSSM lightest GP-even Higgs boson 
h in the decoupling limit, where it has almost exactly the properties of the standard Higgs 
particle. The contributions of the W and top quark loops to the hZy coupling are the 
same as in the SM, but additional contributions are induced by chargino and top squark 
loops. In the low tan/? scenario where the h boson is lighter than ~ 80 GeV, these 
contributions can induce large deviations of the Z ^ h'y decay width from the SM value. 
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even in the decoupling limit. This is illustrated in Fig. 7, where the deviations due to 
charginos and i quarks are shown for two values of the masses 100 and 250 GeV. As can 
be seen, in very large areas of the MSSM parameter space, the deviations can exceed the 
level of several percent. 

If future linear colliders could spend a few months running at the Z resonance, 
a large sample oi Z ^ 7 +Higgs decays could be collected with the expected high- 
luminosities, if the decay is kinematically allowed. In this case, the Higgs-Zy coupling 
could be measured with a precision at the percent level, allowing a stringent test of the 
coupling. In the Standard Model, the couplings of the Higgs particle to W bosons and 
top quarks can be measured with a good precision. In the MSSM, since the contributions 
of the genuine SUSY particles to the decay width exceed the percent level in large areas 
of the parameter space, the h boson can be distinguished from the standard Higgs boson 
even in the decoupling limit. The measurement of the Higgs-Zy coupling in the decay 
Z y+ Higgs will be in this sense, competitive with the measurement of the Higgs 
couplings to two photons at high-energy yy colliders. 
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Figure Captions 

Fig. 1: Feynman diagrams contributing to the Higgs coupling to a photon and a Z boson in 
the SM (a), additional contributions in the THDM (b) and additional contributions 
in the MSSM(c). 

Fig. 2 : The W and top quark form factors in the SM as a function of the Higgs boson mass 

(a) and the branching ratios of the decays hf® —Zj and Z in the 

SM. 

Fig. 3 : The amplitudes for the contribution of the W boson loop and of the sum of the t, b 
loops as a function of (a) and the contribution of the charged Higgs boson loops 
as a function of (b). The contributions are in a THDM with tan/? = 1.6,5 
and 50. 

Fig. 4 : The amplitudes for the contribution of the slepton (a) and squark (except stop) loops 

(b) as functions of the loop masses for tan/? = 1.6,5 and 50. We have neglected 
sfermion mixing and have set Ma = 1 TeV. 

Fig. 5 : Contours in the plane, for which the contribution of the top squark loops 

to the hZ'j coupling is \Ai\ = 0.5, 0.2 and 0.1 for Ma = 1 TeV (a) and Ma = 80 
GeV (b) with tan/? = 1.6. 

Fig. 6 : Contours in the {M 2 , fi) plane for tan/? = 1.6 and Ma = 1 TeV (a) and Ma = 80 
GeV (b) for which the contribution of the chargino loops to the hZ'y coupling is 
\A^\ = 0.2, 0.5 ,1 and 2. Also included are the contours for which the lightest 
chargino mass is m + =70. 

Xi 

Fig. 7: The deviations of the h Z 7 decay width from the Standard Model value [in %] 
for tan/5 = 1.6 and the loop masses = 100 and 250 GeV. (a) Deviations due to 
the chargino loops as a function of M 2 for both signs of /i, and (b) deviations due 
to the top squark loops as a function of 
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